On Axion's Effect on Propagation of Monochromatic Electromagnetic Wave Through 

Strong Magnetic Field 



(N 

o 

(N 
IX, 



i 

-i— > ■ 



> 
00 

o 

O 

t— 1 

o 



Mikhail Khankhasayev and Carol Scarlett 
Physics Department, Florida A&M University 
(Dated: February 8, 2012) 

A possibility of detecting the effect of photon-axion mixing in a cavity experiment is discussed. 
There are two photon-axion modes that acquire different indices of refraction and split in an inho- 
mogeneous magnetic field. For a magnetic field inhomogeneous in the direction transverse to the 
light propagation an analytical solution is obtained both for the index of refraction and the beams' 
trajectories. In a cavity experiment, the beam splitting creates a bifurcation effect, which results in 
a decrease of the light intensity in the central region. Modulation of magnetic field can separate this 
effect from background by providing a narrow frequency range for any observed signal. When one 
integrates this effect over time and accounts for bandwidth, the overall drop in FWHM intensity is 
of order 10 _2 %. This is a very measurable effect. 



I. INTRODUCTION 

It is well known that magnetic fields change optical 
properties of the vacuum (permiability and pemittivity) 
due to photon polarization effects (see, e.g., [1|). The 
presence of axions creates additional optical effects via 
mixing of photons and axions in magnetic field. In ad- 
dition to changes of the index of refraction the axion- 
photon interaction creates a Stern- Gerlach type effect 
of beam splitting in an inhomogeneous magnetic field 
(see, e.g. @, Q. The possible existence of axions have 
been the source of much theoretical and experimental in- 
terest. The axion was first introduced in Refs. [1, [B[ 
to solve the strong CP problem. In addition, as a light, 
weakly interacting and long-lived particle, the axion is 
considered as one the top candidates for explaining the 
dark matter (CDM) in the Universe (see, e.g., Refs. 0- 

i- 

In Rcf. 3], the possibility of detecting the mixing of 
axions and photons by studying the arrival times of ra- 
dio pulses from highly magnetized pulsars was discussed. 
In the present letter we discuss the possibility of detect- 
ing the photon-axion mixing effects in a cavity experi- 
ment. The paper is organized as follows. In Section 2 wc 
present a theoretical framework describing propagation 
of monochromatic light through a magnetized region. In 
Section 3 we present a discussion and quantitative esti- 
mates of measurable effects induced by axions. The final 
section summarizes our results and conclusions. 



II. THEORY 

The description of the photon-axion interactions in an 
external magnetic field B e is described by the Lagrangian 



density 
L = - 



\F^ V + \{d^a - m a 2 a 2 ) + ±g a aF^F^ 
, 2 



90m e 



(1) 



where F M „ is the electromagnetic field tensor, F^ v = 
\^iiv&~ i F llv is its dual tensor, a and m a are the axion 
pseudoscalar field and its mass, respectively. The third 
term describes the photon-axion interaction, where g a is 
the axion-photon coupling constant. The fourth term 
describes the photon's dispersive effects induced by the 
external magnetic field [![. 

In the presence of an external magnetic field it is con- 
venient to present, following [l|,the magnetic field as the 
sum, B — > B + B e .The photon-axion interaction term 
due to the external magnetic field is given by 



Lint = \g a aF^F^ = g a a(B e E) 



(2) 
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In the present paper we limit our considerations with a 
static (or slowly varying in time) magnetic field directed 
along the x-axis, i.e., B e — e.\B e — (B e ,0,0), and uni- 
form along this direction. However, it is inhomogeneous 
in transverse direction (chosen as the y-axis) to the pho- 
ton beam propagation. This conditions correspond to a 
typical setup in a cavity experiment when a laser beam 
propagates perpendicular to the magnetic field. 

According to Eq.@, only the x-component of the pho- 
ton's electric field contributes to the photon-axion inter- 
action, and, for simplicity, we will consider linear, polar- 
ized light along the magnetic field, i.e., E — e\E\ — 
(Ei,0, 0) = —dAi/dt, where A\ is the corresponding 
component of the vector potential. In addition, we as- 
sume that the mirrors' surfaces are parallel to the ex- 
ternal magnetic field. Therefore, during the reflections 
from mirror-to-mirror there will no change of the initial 
polarization of the light. Finally, the initial direction of 
the light beam is chosen to be along the z-axis. Since the 
photon-axion interaction is very weak, this will be the 
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dominant direction for the propagation of light through- 
out the experiment. 

To derive the system of equation describing axion- 
photon interaction in an external magnetic field let's ne- 
glect for a moment the photon-photon interaction term 
in Eq.((T]). The Lagrange equations in this case are the 
following: 



^1- 



(□ - m a 2 )a - g a B, 



dAi 
dt 



= 0, 



(3) 



where □ = V 2 

To study propagation of photons and axions through 
the magnetic field it is naturally to look for the solutions 
in an eikonal form, 



A 1 = Ae^-"^, 



(4) 



where A and a denote the amplitudes of the photon and 
the axion correspondingly. By using this Anzatz we as- 
sume that both the photon and the axion beams prop- 
agate in space and pass through the external magnetic 
filed region together. 

Substituting these expressions in© we obtain, 

(uj 2 — k 2 )A — iujg a B e a = 0, 

(uj 2 - k 2 - m a 2 )a + iuig a B e A = 0. (5) 

Here, we would like to note that in Refs.(@,i] the imag- 
inary unit in the interaction term is missing. 

As it was shown in the photon-photon interaction 
in the presence of an external magnetic field creates dis- 
persion effects which can be taken into account by intro- 
ducing an effective photon's mass Q 7 in Eq.([5l 



(uj 2 - k 2 + Q n )A - iuig a B e a = 0, 
(uj 2 - k 2 + Q a )a + iujg a B e A = 0. 



Here, Q a 



and 



7a , B e 



Q 7 = w 2 -( 



■) 2 , 



(6) 



(7) 



457r B crit 

where B cr u = m 2 /e sa 4.4x 10 13 G is the critical magnetic 
field strength [2j. It is important to note that Q 1 is a 
positive quantity. 

This system of equations can be written in a matrix 
form as 



(L0 2 -k 2 )I + 



™2 



= 0, (8) 



where Qm = ug a B e , and / is the unit 2 x 2-matrix. 
Introducing A = k 2 — u> 2 , we can present this system 

as 



iQm 



A -IQm 
Qa-X 





'A' 




a 



= 0. 



(9) 



The eigenstates and eigenfuctions of the photon-axion 
system in the external magnetic field can be found by 
diagonalizing a Hermitian matrix 

Qy — A —iQm 

iQM Qa - A 

The eigenvalues are give by 



X± = \{Ql+Qa)±\l\{Q 1 



?a) 2 +Q 2 M , (10) 



or 



k 2 ± 



r -(Q-y + Qa) ± \l-AQy - Qa? + Q^ll) 



The matrix R that mixes photons and axions in the mag- 
netized area is given by 



coscp ism<j 
isind) cosd 



where angle <p is determined by 
-tan20 = — - 



M 



(12) 



(13) 



The photon-axion eigenstates A' and a! are given by 

A 1 = coscpA + is'm(f>d 

a' = coscfra + ismcj)A (14) 

in terms of pure photon A and axion a states. In the vac- 
uum without magnetic field (let's label it below as Re- 
gion l)these amplitudes satisfy to the free Klein-Gordon 
equation. 

A. Birefringence 

A physical sense of the obtained eigenfunctions can be 
gained by considering the limit of Qm — > when the 
magnetic field is zero. In this limit, A' — > 
tildeA and a' — > a. It means that the photon and axion 
entering a magnetized region (Region 2) will form co- 
herent mixed photon-axion states propagating with the 
indices of refraction (see, Eq. (fTTj) . 



nl = 1 



± 



i 

1 

J 2 



Qa) 



Q 2 



(15) 



where " -I- " stands for the 7' mode and " — " for the a' 
mode. 

There is a special case when the mixing of photons and 
axions reaches it's maximum. It happens if in Eq. (|13j) . 



Q 



Qy Qa 



OO, 



(16) 



which means that Qm » Q-y — Q a, or Q 1 = Q a . How 
ever, taking into account that Q 7 > and Q a 



m 2 < 



3 



0, the condition of Q 7 = Q a can never be met in the 
vacuum (see, also the discussion in Ref.Q- 

The mixing angle in the maximum mixing case is <fi = 
45 deg, and the eigenstates are given by 



1 



A' = —{A + i~a) 
a' = -j=(a + iA). 



(17) 



The indices of refraction for the symmetric case are given 
by 



2 i , Qm 

n± = 1± — 5-, 



(18) 



(19) 

where /3 = g a B e /2uj. In the symmetric case the indices 
of refraction are linear in the magnetic filed strength. 

The symmetry between the photons and axions in the 
magnetic field was discussed in [9] in the context of the 
duality of photons and axions in the zero axion's mass 
limit. The symmetry is exact under the following condi- 
tions only: (a) the axion's mass is zero, and (b) the pho- 
ton polarization effects in a magnetic field are ignored. 
Otherwise, the duality between photons and axions is 
approximate. 



B. Beam splitting effect 

Calculated above the indices of refraction, Eqs. (fT5| . 
give the light propagation trajectories through the mag- 
netized region using well-known equation (see, e.g., Ucj), 



dl 
Jl 



Vn - l(lWn) 



(20) 



where I is the unit tangent vector along the light propa- 
gation path, and dl is the infinitesimal element of length 
along the light line. From this equation it follows that in 
the area described by a constant index of refraction corre- 
sponding to a uniform magnetic field the light will propa- 
gates along a straight line. Since the indices of refraction 
are different for the photon induced mode and the axion 
induced mode in accordance with Eq. (|15[) the light beam 
entering from Region 1 (free space) into Region 2 (mag- 
netized area) will be separated into two beams. This 
effect disappears when the incident ray is normal to the 
surface separating Regions 1 and 2. 

In an inhomogeneous magnetic field the indices of re- 
fraction (Ti"5j) are not constant as well. In this case the 
light ray splits into two beams (even if the incident ray is 
normal to the surface). Let's assume, for simplicity, that 
the external magnetic field B e is inhomogeneous in the 
transverse direction which is chosen as the y-axis. In this 
case the index of refraction in Eg. (1201) will depend on y 



variable only, and the system of equations is reduced to 
the equation 



dl n dy v 



(21) 



since the z-component of the unit tangent vector I can 
be found from the condition, l y + l\ = 1. 

Represent the beam's trajectory as z = z(y) it is easy 
to find the following exact solution: 



y/1 - C/n 2 , n > 1, 



(22) 



Vl - Cn 2 , n < 1. 



Here, C is an arbitrary constant to be determined by the 
initial condition l y = l y (y ,za). For example, if the light 
ray enters Region 2 at point P = (yo,Za) normally, i.e., 



l y = 0, then 



1 - r$ n>n , 



-\ 1 - ri n<n , 



(23) 



where n = n(y , z ). 

The trajectories of photon-axion modes in the magne- 
tized region are determined by the following integrals: 

f y n dy 
z(y) = / 7 9/ - y + z ,n> n ; 
Jyo V n \V) - n o 

<V) = -I V / n ^ d l+ Z0 ,n<n . (24) 



yo y/ n o- n2 (y) 



Equations (|23|) and (|24)) present the complete solution 
to the problem of determining both the angle of refrac- 
tion and the path of the light beam passing through the 
magnetized area. 

If the index of refraction has a linear dependence on y, 
i.e., 



n{y) = n + b(y - yo), 



(25) 



where b = dn/dy and tiq — n are taken at y = yo; the 
integrals ([24]) can be calculated analytically: 



z{y) = z + — In — - — ■ r, 

o no{n + n ) 



n > no, 



and 



z(y) = z - \ \j n l- " 2 > n < n o- 



(26) 



(27) 



III. QUANTITATIVE ESTIMATES OF 
MEASURABLE EFFECTS INDUCED BY 
AXIONS 

In this section we limit our analysis with the light beam 
entering Region 2 normal to the surface dividing Regions 
1 and 2 at the point y = 0, z — 0. To separate the 
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effect caused by the inhomogeneity of the magnetic field 
it is convenient to present the index of refraction in the 
following form 



no ± dn, 



(28) 



where uq = n(0, 0) is the constant piece of the index of 
refraction; Sn is its variable part; " + " corresponds to 
the 7' mode (n > 1) ,and "-" to the a' mode (n < 1), 
respectively. One can see that the constant part no is 
different for 7' and a' modes. 

Due to an extremely weak coupling ofthe photon to 
the axion, g a ~ 10~ 10 — 10~ 14 GeU _1 , the axion induced 
effects are extremely small, i.e., Sn -C n . Using Eq. (|2"3")l 
one can obtain the following lowest order approximation 
for l r : 



l y s» ±^/25n/no. 



(29) 

For a linear inhomogeneity (|25j) this equation reads as 
l v « ±^/2by/n . (30) 

For the same reason the shift of the light beam in the 
y-direction, y — l y L is much smaller the corresponding 
distance (L) traveled by light in the ^-direction. Using 
Eqs. (j2"o1 |2"7| one can find the following relationship be- 
tween these two quantities 



2n 



(31) 



Substituting this equation into (|29p we obtain 

l y = sin6 =>■ 6 = ±bL/n Q , (32) 

in terms of the angle of refraction 9. The splitting angle 
Ad = + — 6- between the photon and axion induced 
modes is then given by 



AO = 2bL. 



(33) 



This formula is derived under the assumption that the 
magnetic field is inhomogeneous linearly in the transverse 
(y) direction, i.e., 

where Bf = dB e /dy\ y=Q . 

For the case of maximum mixing of photons and axions 
determined by Eq. (|3"4"|) we obtain: 



n ± rj l±/3 , 



(34) 



where /3 = g a BQ/2uj, is the constant component of 
the external magnetic field, giving 



9aB\L 
2wno 



(35) 



Therefore, the splitting angle AO = 6 + — 9_ between two 
modes is 



AO 



9aBfL 



ui 



(36) 



The obtained formula similar to that of Ref. 
9afcB e /uj, where fa is so-called geometric factor. Com- 
paring these formulas, one obtains 

ri In R e 

fa = {B\j B e )L = —— L, (37) 



dy 



where L is the distance travel by the photon beam along 
the z-axis in the magnetic field with gradient Bf along 
the y-axis. 

In Ref. [H, the authors discussed the possibility to de- 
tecting axions by studying radio pulses arriving from 
magnetars with magnetic fields of the order of 10 16 G and 
the radii ~ 10 5 m. To estimate the effect the authors used 
fa = 0.1, which translates to the magnetic field gradient 
of Bf ~ 10 11 G/m. This strong variation of the magnetic 
field over distances comparable to radio wave lengths can 
create an observable splitting angles ~ I0~ 2 rad. In labo- 
ratory experiments the typical strength of magnetic field 
is ~ 10 4 — 10 5 G, the size of the magnetic field region 
is ~ lm, and the magnetic field gradient is ~ 10 6 G/m. 
The geometric factor in this case is fa ~ 50. Assuming 
a photon's wave length ~ 10 _6 m and using (|36[) we can 
obtain the the following expression for the splitting angle 
AO ~ 10 _5 <? a as a function of g a measured in GeV^ 1 . If 
we take the CAST limit [llj for the upper value of g a , 
which is 10~ 10 then we predict a splitting of ~ I0~ 15 rads. 
One way to build up the effect, see Eq. (|3"6")l . is to increase 
the length L of the field region. In a cavity experiment, 
allowing the light to bounce back and forth between two 
mirrors will effectively increase the magnetic field length 
up to about 10 5 m, where mirror reflectivity reaches its 
limit. Multiplying our effect by this total distance would 
give a splitting angle of ~ 10~ 10 rad. HOWEVER, for a 
cavity experiment, reflection from each mirror destroys 
the axion-photon modes. After the first bounce, two 
beams reenter the magnetic field and divide giving four 
new beams, then eight, etc. Importantly, the polariza- 
tion of each new beam remains unchanged. The effect is 
a bifurcation at the face of each mirror and in a distance 
of only 10 4 m, less then mirror limitations, results in an 
average separation of the beam intensity of ~ I0~ 9 rad. 
Assume for simplicity that two planar mirrors enclose a 
region of the magnetic field gradient that is ~ lrn in 
length. Figure QJhows a sketch of a bifurcating distribu- 




FIG. 1. Bifurcation diagram: thick black lines represent the 
bifurcating rays, two bold lines represent a linear split of a 
beam 

tion for only a few bounces. The two thick, straight lines 
show the linear solution for comparison. The splitting 
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angle has been exaggerated to illustrate how much more 
rapidly a bifurcation diverges than a linear separation. 
It can be shown that the overall density of rays for the 
bifurcated distribution expands at a rate of Az^ . Fig- 
ure [5] uses this formula to show that after 10 4 bounces 
assuming the above typical cavity experimental parame- 
ters, (9 ~ 10~15), the rays are shifted away from center 
by as much as 10 _9 m. In Figure [2j Weighted Position 
refers to the position of the new rays weighted by the frac- 
tion of initial light intensity each represents. Summing 
the weighted position gives the average separation of the 
light intensity from the central region. For a Gaussian 
beam, energy spreads away from the center and is, for 
such a small shift, roughly linear giving rise to a drop in 
the FWHM of - lCT 9 ^, where E is the initial energy. 
Furthermore, we can suppress backgrounds through mod- 
ulation of the external field. When one integrates this 
effect over time and accounts for bandwidth, the over- 
all drop in FWHM intensity is of order 10~ 4 , which is a 
very measurable effect. A more detailed computer simu- 
lation results and computer code for the analyses of the 
bifurcation process that incorporates both the realistic 
magnetic field configuration (e.g., for a quadrupole mag- 
net), shape and position of mirrors, etc. will be presented 
in the subsequent paper. 



1E-5-= 

ir-6 n 



S 1E-8 n 

o 

Q. 

"S 1E-9- 

g 1E-10, 



IV. CONCLUSION 

The present paper was focused on exploring a possi- 
bility of observing the mixing of axion and photons in a 
cavity experiment with an inhomogeneous magnetic field. 
The discussion was heavily based on the results of Refs. 
( @> 9 ) • The presented formalism provides a consistent 
way of calculations the effects related to the mixing of 
photons and axions passing through the magnetized area. 
Eq. ([5]) and the obtained mixing matrix (fT2)l correctly al- 
low one to analyze the photon-axion mixing effects at 
arbitrary mixing angle providing a correct mixing phase 
and continuous transition to the symmetric mixing case 
corresponding ip = 45°. For an external magnetic field in- 
homogeneous in the direction transverse to the incident 
light we obtained an analytical solution, Eqs. (PZ3)l and 
(f!M|). for the index of refraction and the beam's path. 
Assuming a linear inhomogeneity of the index of refrac- 
tion we obtained expressions, Eg. and (PZT|) . for the 
trajectories of the photon and axion induced modes in a 
magnetized area. The formulas provide an effective tool 
for analyzing both qualitatively and quantitatively prop- 
agation of light in a cavity experiment since they can be 
used to calculate both the direction and the point of entry 
into the magnetic field. Axions in the presence of exter- 
nal magnetic fields change the optical properties of the 
vacuum in several ways. First, the axion-photon mixing 
creates two indices of refraction, Eq. (fl~5)) for the photon- 
and axion-induced modes which creates the effect of bire- 
fringence. In addition, in an inhomogeneous field these 
two mixed photon-axion modes are subjected splitting - 
a Stern- Gerlach type effect. The splitting effect creates 
a bifurcation effect for the light beam bouncing between 
two (or more) mirrors in a cavity type experiment. It is 
shown that bifurcation strongly increases the signature 
of the axion induced effects. 



1 100 10000 1000000 

Distance (m) 

FIG. 2. The spreading of rays for a bifurcating distribution 
(solid line) compared to that of a linear propagation of two 
beams (dashed line). 
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